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Abstract 

> 

■ Recent results from SuperKamiokande and SNO experiments have set 

severe constraints on possible mixings of a light sterile neutrino, u s , with 
the three active species required for a simultaneous explanation of the solar, 
atmospheric and LSND neutrino oscillation data. A consistent scheme has 
emerged from a global analysis of the data wherein two of the neutrinos are 
nearly degenerate with a mass of order 1 eV, that mix significantly with the 
^h. two lighter states. We present realizations of such a mixed 2+2 oscillation 

scenario based on L e + — L T — L s symmetry (Lj stands for the ith lepton 
number) . Breaking of of this lepton number symmetry by a small mass term 
for v s leads to the required large mixings for both the atmospheric and the 
5_i ■ solar neutrino oscillations. Sum rules for the neutrino oscillation parameters 



are derived within this scheme, and are shown to be consistent with present 
data. These models predict U e 3 ~ 0.02 — 0.03, which can serve as a test 
of this idea. We also present gauge models based on mirror extensions of 
the Standard Model that naturally lead to a light sterile neutrino with the 
required mixing pattern. 



I. INTRODUCTION 



Neutrino oscillation data from three different classes of experiments, viz., atmospheric [|TJ, 
solar and LSND J7|, provide conclusive evidence in favor of non-zero neutrino masses 

and mixings and thus for physics beyond the Standard Model. In order to understand 
the specific nature of new physics, one has to determine the detailed pattern of mixings 
and masses that fit various observations. As more and more experimental data become 
available, these details are getting less fuzzy || and perhaps more interestingly, some simple 
possibilities are being either heavily disfavored or eliminated. In this note, we study the 
nature of new physics and possible new symmetries implied by the three pieces of observation 
taken simultaneously. 

Broadly speaking, since the three different types of experiments are sensitive to three 
different scales of oscillation lengths, or equivalently three different values of mass splittings 
Am 2 , a simultaneous explanation of all data requires that we go beyond the conventional 
framework of three known light neutrinos (z/ e , v^, v T ) and postulate the existence of a fourth 
light sterile neutrino, v s . This is the four neutrino picture that we will seek to understand 
in this paper. Needless to say that the need for the sterile neutrino will be somewhat 
tentative until the LSND results are confirmed by new experiments such as the MiniBOONE 
experiment at Fermilab. 

There are three classes of models that can potentially explain all the data: (i) 2+2 
scenario, (ii) 3+1 scenario and (iii) the mixed 2+2 scenario. 

2+2 scenario: 

In this scenario P,p"0| , two of the neutrino are nearly degenerate and have masses approx- 



imately equal to \j Am LSND ~ 1 eV, while the remaining two states are much lighter. The 
two heavy eigenstates are generally chosen to consist mainly of the near-maximally mixed 
weak eigenstates and v T with their mass difference squared of order Am 2 tm ~ 3 x 10 -3 
eV 2 . Clearly this is designed to explain the atmospheric neutrino data via — v T oscilla- 
tion. The light eigenstates are assumed to contain mostly v e and v s with a mass difference 
squared equal to the solar splitting Am 2 . There is a small mixing between the two pairs 
which provides an explanation of the LSND results. 

3+1 scenario: 



In this scheme fill, the three active neutrinos v e ,a,T have masses < J Am„ ta ~ 0.05 eV, 



while v s is heavier with a mass of order J Am 2 LSND ~ 1 eV. The solar and atmospheric 
neutrino oscillations then involve only the active neutrinos and LSND data is explained by 



indirect oscillations involving the sterile neutrino ||12j| . The goodness of fit to the LSND 
data, while acceptable, is not great Only certain values of the Am 2 LSND are allowed 

and that too at the 95% confidence level in this scheme. 

Mixed 2+2 scenario: 

This scenario is a variation of the 2+2 scheme which is necessitated by recent solar 
neutrino data from SNO experiment in conjunction with SuperKamiokande data. Prior to 
the SNO results, the 2+2 scenario where solar neutrino oscillations were assumed to involve 
only the v e and v s was considered to give a good fit to the data although the fit to the energy 
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distribution was not as good as the fit to the rates. (The 3+1 scenario gave a better fit to 
the energy spectrum as well as the rates but could fit LSND data only at the 95% CL.) 

The SNO results have had a major impact on these conclusions. It showed that there 
is a clear gap between the neutrino flux measured by the SuperKamiokande experiment and 
the charged current flux measured by the SNO experiment. If the solar neutrinos oscillate 
only to the sterile neutrinos as is assumed in the original 2+2 scenario, the fluxes measured 
by SNO and SuperKamiokande should be equal. The observed gap in the flux would appear 
to exclude the 2+2 scenario. The mixed 2+2 scenario circumvents this problem. In this 
scenario, solar v e is assumed to oscillate into a linear combination of v s and v^ T . This 
would allow an understanding of the gap in the flux measured at SuperKamiokande and 
SNO, since the v^ T components in solar neutrinos will contribute to neutral current cross 
section at SuperKamiokande, but will not play any role in the charged current cross section 
at SNO. The challenge then is to achieve such a mixed oscillation scenario without upsetting 
the other constraints from accelerator data on the one hand and the atmospheric neutrino 
oscillation data on the other. Atmospheric neutrino data do set severe constraints on such 
a mixed scheme since that prefers — v T oscillations with only a limited amount of v s 
component allowed. 

In a recent paper, Gonzales- Garcia et. al. |T3[ have studied this issue in quantitative 



detail. They proposed a 4 x 4 neutrino mixing matrix where U e s ~ ?7 e 4 ~ was assumed 
in order to satisfy bounds from reactor neutrino data, especially from the CHOOZ [|TJ] and 
PALO-VERDE [[UJ experiments. Here the two heavy eigenstates with nearly degenerate 



mass of order 1 eV are denoted by z/ 3 and while v\ and v 2 denote the lighter eigenstates. 
There is substantial mixing between the heavy and the light sectors. Solar neutrinos can 
oscillate into u^ T with a probability proportional to (1 — \U s i\ 2 — |t/ S 2| 2 ) where U s \ and U s2 
parametrize the non-negligible sterile neutrino component in the two light eigenstates v\ and 
v 2 . This helps reconcile the SNO and the SuperKamiokande solar neutrino data. On the 
other hand, this mixing pattern allows atmospheric z/ M 's to oscillate into sterile neutrinos 
with a probability of 2(|f/ s i| 2 + |f/ S 2| 2 )(l — |^3| 3 ~~ I ^7*4 1 2 ) • The limit on the sterile neutrino 
fraction in the atmospheric data therefore tends to reduce the amount of the needed active 
neutrino fraction in the solar neutrino data. A delicate balance is needed and the authors 



of Ref. [14] found a fit to both the solar (including SNO data) and atmospheric data for 



a range of values of the mixing angles consistent with all other observations. This mixing 
pattern is an interesting way to accommodate the LSND data for a wider range of Am\ SND 
than the 3+1 scenario and in fact this delicate balance means that this mixing pattern can 
be tested once the data improves and as new experiments give results. 

In this paper, we explore theoretical scenarios that can lead to mixed 2+2 scenarios 
with mixing patterns close to the one just described. We find that this can happen if there 
is an approximate leptonic symmetry L e + L M — L T — L s in the Majorana mass matrix 
for the four neutrinos in the context of a seesaw mechanism for neutrino masses. If this 
symmetry is broken by a mass term for u s , the large mixing angles needed for both solar and 
atmospheric neutrino oscillations will result. We show that this scenario predicts certain 
sum rules relating neutrino oscillation parameters which are in agreement with the data 
currently. We then show that if the Standard Model is duplicated to have a mirror sector, 
the lightest mirror neutrino can play the role of the light sterile neutrino with the required 
mixing properties. The desired leptonic symmetry emerges at low energies, if the original 
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model at high energies has the symmetry (L e + L M — L T ) x (L' e + L'^ — L' T ) (here prime denotes 
the symmetries of the mirror sector). 

We have organized this paper as follows: In section II, we describe the phenomenology of 
the mass matrix with the approximate leptonic symmetry described above; in section III, we 
study the phenomenological implications and numerical fits to data for this mass matrix and 
derive certain sum rules for oscillation parameters. In section IV, we present a gauge model 
based on the existence of a mirror sector that leads naturally to a light sterile neutrino with 
the desired mixing pattern. In section V we conclude. 



II. APPROXIMATE L E + L„ - L T - L s SYMMETRY AND THE MIXED (2+2) 

NEUTRINO OSCILLATION SCHEME 



Consider the four-neutrino system (v e , u T , u s ). Let the mass matrix be of the form: 
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(1) 



This matrix has an approximate U(l) symmetry which can be identified as L e + L^ — L T — L s 
where Lj stands for the ith lepton number. This U(l) symmetry is not exact, it is broken 
by the v s mass term {6m§){y s v ^) in Eq. (1). This breaking will turn out to be small, in 
fact 5 will be the smallest entry in Eq. (1). In the limit 5 — > 0, M v of Eq. (1) will collapse 
effectively to a 2 x 2 matrix, with the four neutrino states forming two Dirac fermions. The 
mass matrix in the «5 — > is given by 



M° 





e 2 \ 
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1 ) mo 
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V a 





(2) 



Note that due to the L e + L M — L T — L s symmetry, mixing occurs only between (z/ e , z/ M ) 
and (u T , u s ) in this limit. We shall be interested in the case where a ~ 1 ^> ei, e 2 . The two 
nonzero mixing angles are then 



tan 9 ST ~ a 



1 + a 2 

The two Dirac neutrinos will have masses given by 



(3) 



m h ~ Vl + a 2 m , 



mi 



m . 



(4) 



'1 + a 2 )m 



o- 



The mass splitting relevant for the LSND experiment is Am 2 LSND ~ m 2 — mf 
We shall choose m h ~ 1 eV so as to explain the LSND experiment, along with ei i2 — 
(2 — 3) x 10~ 2 to fulfill the LSND mixing angle requirement. At this stage, no mass splitting 
other than that for LSND is induced. We can allow for the possibility that either e\ or e 2 is 
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zero due to some flavor symmetry. These special cases will reduce the number of parameters 
by one. We shall keep both terms to be nonzero to be general, but in Sec. Ill, we shall also 
discuss these two special cases. 

Now, let us include the effects of L e + L^ — L T — L s symmetry breaking through the mass 
term 5 (the (4,4) entry of Eq. (1)). This mass term, which breaks the symmetry rather 
economically, will serve several purposes. It will induce two more mass splittings, to be 
identified with the atmospheric mass splitting (Am^ tm ) and the solar mass splitting (Atoq). 
At the same time the 5 term will lead to maximal or near maximal mixing between the two 
would be Dirac states both in the atmospheric neutrino sector and in the solar neutrino 
sector. 

On the theoretical side, we envision that the global L e + L^ — L T — L s symmetry is broken 
by some high scale physics. For example, quantum gravity is suspected to break all global 
symmetries, so that could be the source of the 5 term. It is quite natural to assume that 
the symmetry breaking effects show up first in the v s mass term, which is a complete singlet 
of the Standard Model. In fact, in the explicit gauge models that we have constructed (see 
Sec. IV), 5 arises through quantum gravity, and the analogous symmetry breaking effects 
are negligible in all other entries of Eq. (1). 

Including the 5 term, the eigenvalues of M v of Eq. (I) can be computed in the ap- 
proximation 1 ~ a ^> e 12 ~ S. Neglecting quadratic terms in e 12 and 5, these masses 
are: 



for convenience. We have tacitly assumed all parameters of M v to be real for simplicity. 
For 5 positive, we have m\ < m 2 < m 3 < m 4 . If 5 is negative, we can rearrange the labels, 
mi *-* m 2 an d itl 3 <-> m 4 , so that the hierarchy m\ < m 2 < m 3 < m 4 is maintained. 
From Eq. (5) the three relevant mass splittings are found to be 




(5) 



Here we have defined 




(6) 



Am 2 r 





ml ~ (2SX) ml . 



(7) 



The leptonic mixing matrix U is given by (to linear order in 5 and e 12 ) 



5 



V2[(ae 1 +€ 2 )(ei -ae 2 ) + 



2(l + a 2 ) 



(5-X)] 




v"x(x+a) 

V2[(a ei +e 2 )( ei -a e2 )+ 2 

'i+a 2 V(x(x+a) 





Here the entries in the first row are U e i, i = 1 — 4, the ones in the second row are U^, the 
third row entries are U T i and the last row entries are U S i. Here we have used the definition 
v a = Y^t=i U a k^k, where a = (e, fx, r, s) denote the four flavors and k = 1 — 4 denote the 
mass eigenstates. 



III. NUMERICAL FITS 

With the mass splittings and the mixing matrix entries in hand, we can now confront 
the model with oscillation data from LSND, solar neutrino and atmospheric neutrino exper- 
iments. 



A. LSND experiment 

From the expression for Am 2 LSND given in Eq. (7), we see that (1 + a 2 )ml ~ (0.2 — 6) 
eV 2 in order to explain the positive results seen by the LSND collaboration. The (z/ M — v e ) 
oscillation probability relevant for LSND is given by 

P^ALSND) ~ 4 \U: 3 U, 3 + U: 4 U„\ 2 sin 2 (^|f^) , (9) 

where we have ignored the contribution proportional to the smaller Am^ and Am 2 tm and 
assumed that m 3 ~ m 4 . Comparing Eq. (9) with the matrix elements of Eq. (8), we make 
the identification 

sin 2 26 LS ND » 4 |[/ e * 3 [/ M3 + t/ e * 4 M 2 - 4 

Here sin 2 29lsnd — 3 x 10~ 3 is the mixing parameter usually quoted in the two flavor 
oscillation analysis. The experimental observation can be explained by choosing (aei + 
e 2)/(l + a 2 ) — 0.03. Note that for this analysis, the breaking of L e + L^ — L T — L s symmetry 
is not significant, we could have obtained identical results by performing the two flavor 
oscillation using Eq. (2) in the exact limit of this symmetry. 

Note that the effective mixing parameter U e3 , normally discussed in the three neutrino 
oscillation scenario, is given in our model by \U^l^\ 2 = (\U e s\ 2 + |t/ e 4| 2 ) since m 3 ~ m 4 . Thus 
|C/g3 | ~ 9 lsnd — (0.02 — 0.03). This is a definite prediction of the model that can be used 
as one of its tests. 



ae 1 + e 2 



(1 + a 2 ) 



(10) 
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B. Solar and atmospheric neutrino oscillations 



We shall follow the global analysis of solar and atmospheric neutrino oscillation data 
carried out in Ref. [14 including the recent SNO results. In that paper it was assumed that 
U e3 = U ei = 0. This approximation holds to a high degree of accuracy in our case, since 
U e3 ~ U e4 ~ 9 LS nd/V2 ^ 0.02. Such a small mixing of v e with the heavier mass eigenstates 
is insignificant for the analysis of solar and atmospheric neutrino oscillations. 

Consider first the atmospheric neutrino oscillations. The parameter space is determined 
by {Ami, IM7(IM 2 + IM 2 )> (l^il 2 + W), (|t/ M i| 2 + U^l 2 )}. Here the mass 
splitting is given in Eq. (7), which we shall choose to be Am 2 tm ~ (1 — 6) x 10~ 3 eV 2 . The 
mixing parameter |t/ M 3| 2 /(|f/ M 3| 2 + |f^t4| 2 ) — 1/2 from Eq. (8). This is the leading parameter 
that controls the disappearance of z/ M through oscillations into either v T or v s . Our model 



prediction agrees quite well with the results of the global analysis [0. In particular, our 
model predicts the deviation of this parameter from 1/2 to be extremely small, of order 
1%. The parameter (|?7 s i| 2 + (t^l 2 ) — 1/(1 + a 2 ) parametrizes the projection of the sterile 
neutrino component in the atmospheric neutrino oscillations. When this parameter is equal 
to 1, we have pure — v T oscillations, while if it were zero we have pure — v s oscillations. 
The atmospheric neutrino data prefers oscillations into mostly z/ r , but significant sterile 
component is also allowed. When combined with the solar neutrino analysis, which prefers 
this parameter to be small rather than large, Ref. |14|] obtains (|f/ s i| 2 + |£^2| 2 ) — 0.21 — 0.5 
as preferred JT3J. That fixes the parameter a = (1 — 2). Lastly, the mixing parameter 
(I^Vl 2 + I^msI 2 ) — @lsnd + o 2 ^ 2 / (1 + a 2 ) 3 , as can be seen with a little algebra. Numerically 
this is about 0.001, which is too small to be of significance in atmospheric oscillations. 
Thus, our model corresponds to the restricted case of (jC^il 2 + |t^2| 2 ) = studied in Ref. 



1411711 . This scenario gives a reasonable global fit with a goodness of fit (GOF) quoted to 



be 59% |H |. Actually, within our scheme this GOF will be somewhat better since one of 
the parameter \U / {\U + l^/tl 2 ) is fixed to be 1/2, which is very close to the central 
value of the experiments. 

Thus, a good fit to atmospheric neutrino oscillation within our scheme requires a ~ 
(1 - 2), and 5 = (2 x 10" 3 - 3 x 10~ 2 ), if we make use of Am 2 LSND ~ (0.5 - 1) eV 2 . The 
parameter (aei + e 2 ) lies in the range (0.06 — 0.15), obtained from fitting the LSND mixing 
angle. We see that 5, ei 2 <C 1 while a ~ 1, consistent with our approximations. Note also 
that while 5 is somewhat smaller than ei ; 2, it is not much smaller. In particular, 5 may be 
comparable to (ei — ae 2 ) for a wide range of parameters. This observation will be relevant 
for the solar neutrino mixing angle prediction. 

Turning to solar neutrinos, the parameter space is characterized by (Am^, |?7 e 2| 2 /(l — 
I L/e2 1 2 ) , (|^si| 2 + I f^s2 1 2 ) ) - The mass splitting is given in Eq. (7). The first mixing parameter 
is recognizable as tan 2 9 & that is usually used in the two flavor oscillation analysis for v e 
disappearance. The second mixing parameter specifies the amount of sterile neutrinos in 
solar neutrino oscillations. When (|C/ s i| 2 + |t/ S 2| 2 ) = 0, v e oscillates only to active species, 
when this parameter is 1, f e oscillates into a pure sterile state. Solar neutrino data prefers 
v e oscillating predominantly to an active species, with some sterile admixture allowed. 

From Eq. (8), we see that tan 2 0© ~ (X - 5)/(X + 5). If S < X, this angle will be 
maximal. However, S and X may be comparable, so tan 2 #© can deviate significantly from 
one. Note that the deviation will be to lower values compared to one, a feature that goes 
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well with oscillation data. 



C. Sum rules for oscillation parameters 

Now we show that the model under study leads to an interesting sum rule when the 
oscillation data from LSND, solar, and atmospheric neutrinos are combined. To see this, 
note that (|£/si| 2 + |£/s2| 2 ) — 1/(1 + a 2 ) fixes the value of a. The solar oscillation angle 
tan 2 #0 fixes the value of 5/X, while the ratio Am|/Am^ SJV£) fixes the product 5X. The 
parameters (a, 5, X) are then completely fixed. The sum rule arises by examining the ratio 
Am 2 atm / Am 2 LSND ~ 2a 2 5/(l + a 2 ) 2 , which only depends on the same set of parameters. A 
simple calculation shows the sum rule to be 

\Am 2 } 2 1 

This sum rule is obeyed by current experiments rather well. As an example, let us choose 
(\U sl \ 2 + \U s2 \ 2 ) = 0.25, Am 2 atm = 3 x 1(T 3 eV 2 , Am 2 LSND = 1 eV 2 , tan 2 6 & = 0.5. Eq. (11) 
then predicts Awig = 9.6 x 10 -5 eV 2 . This value is nicely consistent with large angle solar 
neutrino oscillations. 

Consider now the special case e\ = 0, e 2 ^ 0, which may be imposed by some flavor 
symmetry. In this case, there is an additional sum rule, which may be taken to be a relation 
for the LSND mixing angle. (This can be seen by noting that 4(1 + a 2 ) (e\ — ae 2 ) 2 ~ X 2 — 5 2 .) 



Q ~ ^""atm K\ u sl\ T \u s2 \ ) tail *(/ Q 

6 LSND ^ TT- 2 M ITTT9. ITT ma/9 " ( 12 ) 



Am 2 tm (\U 8l \ 


2 + l 


u s2 \ 


2 ) tan 20 


4Am| SiVD (1 - 


\u sl 


2 _ 


\Us2\ 


2\3/2 



To see how this sum rule compares with experiment, let us take Am 2 tm = 2 x 10 3 eV 2 , 
Am 2 LSND = 0.5 eV 2 , (\U sl \ 2 + \U s2 \ 2 ) = 0.7, tan 2 6 & = 0.75. From Eq. (11) this choice will 
predict Am| = 4.4 x 10~ 4 eV 2 , and from Eq. (12), 9lsnd = 0.03. Both predictions are in 
good agreement with observations. 

In the opposite case, where e 2 = 0, e\ ^ 0, there is again a sum rule analogous to Eq. 
(12), given by 



? _ Am lt m tzn29 Q 

LSND - AAm 2 LsND (1 _ ^2 _ |[/ s2 |2)l/2 



(13) 



In this case, Olsnd tends to be small, Olsnd — 0.013 for the same set of parameters as 
above. So this special case appears to be disfavored by present data, unless LSND mixing 
angle settles to a much smaller value. 



IV. GAUGE MODELS 

The mass matrix of Eq. (1) can be obtained from extensions of the Standard Model that 
incorporate the seesaw mechanism along with an approximate L e + — L T — L s symmetry. 
In order to have a naturally light sterile state v s in the light spectrum, we resort to the 
idea of mirror universe |18|JT9|] . In this scenario, it is assumed that there is a mirror sector 



S 



to the Standard Model, which is identical to our world in terms of its particle and force 
content. The symmetry breakings can however be different [[HJ. The main advantage of this 
extension is that the existence of the extra gauge quantum numbers of the new neutrino 
species makes it easier to understand the smallness of their masses. For instance one can 
employ a type II seesaw mechanism, i.e., seesaw via tiny VEV of a superheavy triplet Higgs 
|20| instead of the heavy right handed Majorana neutrinos in both sectors to make all active 
as well as sterile neutrinos light. The mirror neutrinos, being Standard Model singlets, can 
play the role of sterile neutrinos and one has a natural way to understand why a singlet 
fermion, u s , may have such a tiny mass. 

As mentioned before, the scale of electroweak symmetry breaking in the mirror sector 
v' is assumed to be different from that in the familiar sector (v) This is achieved by 
incorporating a mirror odd field r\ which has a VEV and asymmetrizes the SU(2) doublet 
Higgs mass terms in both sectors. ( Equivalent ly, this can be achieved by the soft breaking of 
the mirror symmetry.) Due to different weak scales, the mirror neutrinos will have different 
masses compared to the Standard Model ones. If there is mixing between the neutrinos 
of the two sectors, which is allowed by gauge symmetry, then oscillations between the two 
types of neutrinos will occur. 

Assume that our world respects L e + L ll — L T and the mirror world respects the analogous 
L' e + L'^ — L' T . If these symmetries were exact, a linear combination of familiar neutrinos, 
u e and Up, as well as another one of mirror neutrinos, i.e., u' e and u' , will remain massless. 
We shall identify the second linear combination as v s . Such global symmetries are perhaps 
good only up to quantum gravitational corrections. We assume that they are indeed broken 
by Planck scale effects. Therefore, once the Planck scale breaking terms are included in 
the Lagrangian, both the massless neutrino states would pick up small masses. They give a 
mass to v s of order 5m ~ v' 2 /Mp\. If we choose the mirror weak scale, v' ~ (2 — 4) TeV, 
one finds m Vs ~ (2 — 8) x 1CT 3 eV. There will be Planck suppressed corrections to other 
entries in the neutrino mass matrix as well, but they will not be of much consequence for 
our purpose. For example, the v e v e entry will be corrected by a negligible amount ~ 10 -5 
eV. 

Let us now discuss how the dominant entries of the neutrino mass matrices arise. In the 
familiar sector, the dominant entries of M u , in units of mo, are the 1 and the a entries. The 
1 arises via the seesaw mechanism that we describe now. 

It proves convenient to use the type II seesaw formula [^] for this purpose, where the 
tiny VEV of a B — L = 2 triplet, rather than a heavy Majorana neutrino, gives small mass 
to the neutrino. To implement this mechanism p0| |, we introduce into the Standard Model 



in both sectors, SU(2) triplets A^ and A' L with hypercharges +2. These fields will have the 
following Yukawa couplings to the lepton doublets: 

£yuk = fijiLji^r.ALLj + LfV 2 r.A' L L;.) + h.c. (14) 

By itself, Eq. (14) does not break lepton number, but the A^ and A' L fields will also have 
gauge invariant couplings to the Standard Model Higgs doublet fields that violate lepton 
number. The origin of the naturally small VEVs can be seen from the following effective 
potential, where mirror symmetry is assumed to be softly broken: 

V(<f>, 0', A x , A' L ) = +M 2 (AiA L + AlA' i ) + ( / iA^ + /i'A2W + (15) 
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Here the ... denote higher order polynomial terms that are not relevant for the present 
discussions. The scales in Eq. (15) is chosen to be high seesaw scales in both the sectors. 
The VEVs of the triplet fields are given by 

The value of M 2 / fi is fixed so that the active species has a mass term (connecting and u T ) 
of order 1 eV (the 1 entry). That sets M 2 / /i~3x 10 13 GeV. The mirror seesaw scale M 2 / fjf 
is not identical, we shall choose it to satisfy the cosmological requirement that during big 
bang nucleosynthesis, no more than one sterile neutrino species be in thermal equilibrium 
with the photons. This requires M 2 /// ~ (10 6 — 10 7 ) GeV. The heavy mirror neutrinos will 
then have a degenerate mass of order (10 — 100) GeV. 

The mixing term amo of Eq. (1) is also of order 1 eV. This entry involves the mixing 
of v e with v a and thus will require the breaking of the two L e + — L T symmetries to a 
single i7(l). This is achieved by choosing a Higgs field x, which transforms like a bidoublet, 
i.e., (2, 1; 2', 1') under the two gauge groups and carries charge of (—1, —1) under the two 
leptonic {7(1) x U(l)' symmetries. Its VEV will break the {7(1) x U(l)' to the diagonal {7(1) 
subgroup {(L e + — L T ) — (L' e + L'^ — L' T )}. This is equivalent to L e + — L T — L s since 
two of the mirror states are heavy. 

The x fi e ld has couplings to the lepton doublets of type L e xL' e . The field x a ls° acquires 
a VEV via seesaw type coupling like the triplets, i.e., V D rf' x (j)(f>' + M 2 x X- Since amo ~ 1 
eV, we need (x) ~ iL nr- ~ 1 eV. This can be obtained by choosing v ~ 200 GeV, v' ~ 2 

X 

TeV, fi" ~ M x ~ 10 14 GeV, so that amo ~ 0.4 eV. As for the entries ei^, they will arise in 
a manner similar to the entries 1 and a except that their magnitudes have to be explained 
by some effective Yukawa couplings being order 0.02. As discussed before, to get a nonzero 
value of 5 ~ (2 x 10~ 3 — 3 x 10~ 2 ), we invoke the Planck scale induced {7(1)' violating 
couplings and v ' ~ 4 TeV. This is about 25 times larger than v. So cosmology of mirror 
sector will be identical to what has been studied already pTJ. 

Since this model has two heavy sterile neutrinos (with masses in the multi GeV range), 
one must address their cosmological implications. Their weak interaction has a strength 
of Gpe ~ 10 -3 GV due to the higher mirror symmetry breaking scale. From this, one 
can calculate their decoupling temperature of the heavy mirror neutrinos by using the 
out of equilibrium condition as follows. The mirror neutrino interaction rate is given by 
G 2 7 e 2 M 2 (T'M) 3 / 2 e~T r where T' is the temperature of the mirror sector when the tempera- 
ture of the visible sector is T, and M is the mass of the heavy neutrino. Let us define the 
ratio of these temperatures to be (3 = T'/T. The out of equilibrium condition then gives 
the following equation: 



n^To) MpGpt 2 M 3 ' 1 ' 

Here %d stands for the ratio M/Td where Tp> is the decoupling temperature, xd is deter- 
mined by the following equation: 

G 2 F 6 2 [3 3 M 3 Xr i 2 e-^ ~ . (18) 

Nip 
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This leads to %d ~ 3.6. Using this, we see that for M = 100 GeV, the contribution of mirror 
heavy neutrinos to the energy density at the BBN epoch is — = n ^ T J^ M ~ 0.036. Thus, 

J OJ J 1 Pv n-f(T D )T B BN 

even with both the heavy neutrinos, the total contribution to the energy density during 
nucleosynthesis is negligible. The lighter mirror neutrino will be in equilibrium during 
nucleosynthesis, so the model predicts N„ — 4 for nucleosynthesis. 

The heavy mirror neutrinos are not stable, they decay dominantly into the light mirror 
neutrino and a mirror photon. The decay lifetime can be estimated to be r~ x {y' — > 1^+7') — 
(a/2)m'3[3eG F m' u m 2 T /(32ir 2 m^)} 2 ~ 10" 4 sec. for m' v ~ 100 GeV. This will help eliminate 
most of the heavy neutrinos from the universe. 

The mirror universe hypothesis is one possible framework that explains naturally the 
existence of a light sterile neutrino. Other approaches based on radiative neutrino mass 
generation mechanisms p2| may also provide realizations of our mass matrix, Eq. (1). 



V. CONCLUSIONS 



We have presented a simple ansatz for the four neutrino mass matrix which leads to a 
mixed 2+2 scenario that fits all neutrino oscillation data (solar, atmospheric and LSND). 
The mass matrix follows from L e + — L T — L s symmetry (Lj is the zth lepton number). 
When this symmetry is broken by the small mass of u s , maximal mixings required for 
explaining the solar and atmospheric neutrino data result. We have derived certain sum 
rules for neutrino oscillation parameters within this scheme. These sum rules are consistent 
with present data. The model predicts a small but nonzero value for the effective U e s entry 
of the leptonic mixing matrix, which is approximately equal to the LSND oscillation angle, 
~ (0.02 — 0.03). This prediction can be used to test the model. We have also presented a 
gauge theory realization of the neutrino mass matrix based on the mirror sector hypothesis. 
This specific realization leads to interesting dark matter cosmology that has been studied 
previously. 
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